Einstein boundary conditions for the 3+1 Einstein equations 
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In the 3+1 framework of the Einstein equations for the case of vanishing shift vector and arbi- 
trary lapse, we calculate explicitly the four boundary equations arising from the vanishing of the 
projection of the Einstein tensor along the normal to the boundary surface of the initial-boundary 
value problem. Such conditions take the form of evolution equations along (as opposed to across) 
the boundary for certain components of the extrinsic curvature and for certain space-derivatives of 
the three- metric. We argue that, in general, such boundary conditions do not follow necessarily 
from the evolution equations and the initial data, but need to be imposed on the boundary values of 
the fundamental variables. Using the Einstein-Christoffel formulation, which is strongly hyperbolic, 
we show how three of the boundary equations up to linear combinations should be used to prescribe 
the values of some incoming characteristic fields. Additionally, we show that the fourth one imposes 
conditions on some outgoing fields. 

PACS numbers: 04.25.Dm, 04.20.Ex 



The advent of sensitive gravitational wave detectors 
such as LIGO [l| and LISA [2| motivates the need for 
accurate templates of gravitational wave emission from 
powerful sources such as the collisions of black holes with 
other black holes or with neutron stars [3| . Even though 
the early and late stages of such merging events can be 
modelled analytically by perturbative methods 0, Q , the 
actual merger phase is highly non-linear and can only be 
expected to be modelled accurately by numerical simu- 
lations. Most numerical simulations of dynamical black- 
hole spacetimes break down well before the relevant in- 
formation can be extracted, for reasons that are yet to 
be pinned down. Numerical analysis on the side, there 
arc many fundamental issues in the Einstein equations 
that may be relevant to the stability of any numerical 
implementation of the binary black hole merger, among 
them: intrinsic ill-posedness of the related Cauchy prob- 
lem, constraint violations, poor physical choices of bi- 
nary black hole data and poor choices of boundary con- 
ditions. The solution to long term numerical instabilities 
of binary black hole mergers is likely to involve an or- 
chestrated control of all these factors, in addition to the 
other elements of general numerical analysis that apply to 
time-dependent problems with boundaries. In this work, 
we point out that even before thinking about a numeri- 
cal integration of the Einstein equations in a region with 
timelike boundaries, the possible implications of the Ein- 
stein equations to the boundaries need to be taken into 
consideration, as a matter of principle. To date, most 
numerical simulations with boundaries have largely dis- 
regarded this issue, as have most mathematical studies 
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of the mixed initial-boundary value problem of the Ein- 
stein equations, with exceptions explicitly noted in the 
following. 

Consider the 3+1 formulation of the vacuum Einstein 
equations with vanishing shift vector, in which the metric 
has the form 



ds 2 



-a 2 dt 2 + jijdx z dx J 



(1) 



where a is the lapse function and 7y is the Riemannian 
metric of the spatial three-slices of fixed value of t. The 
Einstein equations G a b = take the ADM form [|| 



-2aK 



Kij = a(R ij -2K il K l j + KK ij )-D i D j a, 
with the constraints 



R 



K^K 10 + K 2 



0, 



D-;K J i - DsK = 0. 



(2) 
(3) 

(4) 
(5) 



Here an overdot denotes a partial derivative with respect 
to the time coordinate d/dt, indices are raised with the 
inverse metric 7 y , Di is the covariant three-derivative 



consistent with 



Ha is the Ricci curvature tensor of 



jij, R its Ricci scalar and K = ^y lJ Kij. Because there 
are no time derivatives of the lapse in any of the equa- 
tions, the lapse function is completely free. Because a 
time derivative of the constraints turns out to be a linear 
combination of the constraints themselves and their space 
derivatives by virtue of the evolution equations, they will 
vanish at subsequent times in the domain of dependence 
of the initial slice if vanishing initially (how stable the 
propagation of the zero values of the constraints is to 
small perturbations does not concern us at this time 

Strictly speaking, the initial-value problem consists of 
finding a spatially periodic or square integrable solution 
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to the Einstein equations in the half-spacetime t > 0. In 
practice, one looks for a spatially non-periodic solution 
in t > 0, x 1 e V, where V is a box in R 3 . There arises the 
problem of specifying data on the timelike boundaries of 
the region of interest. The problem is rooted in the fact 
that in order to specify such data one needs to know the 
values of the solution on the boundary for all time, which, 
of course, one does not know. In an initial value problem 
without constraints, the boundary values may or not be 
constrained by the choice of initial data. It seems to 
us that if the problem, in addition, has constraints, the 
boundary values will know it, and it should be natural 
to expect consequences to the boundary values from the 
presence of the constraints. The point has been raised 
in a precise sense by Stewart Q and further pursued in 
a similar direction by Calabrese et al [j| and by Szilagyi 
and Winicour 10]. 

The existence of at least some such consequences is 
intuitively predictable from the point of view that the 
3+1 equations derive from a 4-d covariant formulation, 
which is more fundamental, that is: G a b = 0. The 3+1 
split is simply a choice of linearly independent combina- 
tions of the ten components of the Einstein tensor. The 
choice is precisely that which removes any time deriva- 
tives of second order from four of the equations, and is 
accomplished by projecting the Einstein tensor along the 
direction normal to the slicing G a i,n b = 0. 

One can easily convince oneself (or else refer to [Til ]) 
that the projection of the Einstein tensor along the di- 
rection normal to a surface given by a fixed value of any 
of the coordinates x % removes all second derivatives with 
respect to that coordinate x l . Let's say we look at our 
boundary at a fixed value of x and write down the pro- 
jection of the Einstein tensor of the metric along the 
normal vector e b = (0,j xl ), which may be transformed 
to unit length by dividing by y/"f xx , but that is clearly 
not necessary. The four components of the projection of 
the Einstein tensor along e b are 

G tb e b = -i 7 " (0n7),« -7 W Tta,i) " KD * a + K k^ k a 



G zb e b 



G x be b — 



a 
K x 



KK X + R X V - -D x D y a = (6b) 



KK X + R X ~ —D x D z a = (6c) 



K K * ~\{R + K^Kn + K 2 ) + KK X 
+R X + - (DW ja - D x D x a) = (6d) 



Here the time derivative of the components of the ex- 
trinsic curvature is applied after raising an index, that 
is: K* = ("/ lk Kk]),t- The reader can verify that R x and 
R x do not involve second derivatives with respect to x of 
any of the variables and that the combination R x — ^R 
doesn't either. In the framework of the ADM equations 



©, even though, to our knowledge, there is no clear no- 
tion of which kind of boundary conditions would be ap- 
propriate for a unique solution, one could think of 16al) 
as an evolution equation for a particular combination of 
the first space-derivatives of the metric, of Eqs. i|6bjl and 



cj) as evolution equations for K x and K x respectively, 



and of Eq. I|6d(l as an evolution equation for the combina- 
tion K*j + K z z . Or one could freely combine them linearly 
to obtain evolution equations for combinations of these 
fundamental variables. In any case, these are evolution 
equations along the boundary itself which can be viewed 
as a sort of mixed Neumann-Dirichlet boundary condi- 
tions [13 

The situation clarifies in the case that one considers 
a first-order reduction of the ADM equations. In that 
case, the space-derivatives of the metric are adopted as 
additional variables Qijk = "fij,k- This has the effect of 
eliminating all appearances of x— derivatives of the fun- 
damental variables from Eqs 10 • However, a straight 
reduction in this fashion is not well posed in the sense 
that it is not strongly hyperbolic. If one is to consider 
a first-order reduction, one is much better off adopting 
one that is strongly hyperbolic in the sense that there 
exists a complete set of null eigenvectors of the principal 
symbol of the resulting system of equations. 

In a generic strongly hyperbolic system of PDE's, 
the null eigenvectors constitute characteristic fields that 
propagate with the characteristic speeds of the system, 
and the notion of consistent boundary conditions can be 
defined unambiguously, as follows [lij. Characteristic 
fields that are outgoing at the boundary are determined 
by the initial values and cannot be prescribed freely at 
the boundary. On the other hand, characteristic fields 
that are incoming at the boundary must be prescribed 
in order for there to be a unique solution, and, further- 
more, can be prescribed arbitrarily (we do not concern 
ourselves with stability for the moment). 

In the case that there exist constraints, it seems to us 
that this notion of consistent boundary conditions may 
still be used but must be adapted to the particularities 
of a constrained system. To wit: if the incoming charac- 
teristic fields are to lead to a solution that satisfies the 
constraints inside the region of interest, then they must 
"know" about the constraints and, in general, may not 
be prescribed as freely as in the case where there are 
no constraints. In the case of the Einstein equations, 
we have shown that there are four boundary equations 
implied by the ten Einstein equations themselves, which 
are completely equivalent to linear combinations of the 
evolution equations and constraints that have no second 
derivatives across the boundary. Moreover, in a sense, 
the boundary equations are to the boundary surface what 
the constraints are to the initial slice. To see this consider 
the Bianchi identities V a G ab = 0, which take the follow- 
ing form in a generic coordinate system (t,x l ) adapted 
to spacelike slices: 



d x G 



xb 



-d t G 



tb 



d y G yb - d z G zb - 



l T a ac G 



cb 



4-rb 



r° ac G a 



(7) 
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where 4 r c a h are the Christoffel symbols of the four-metric 
gab. The right-hand side obviously contains no deriva- 
tives with respect to x higher than second order. There- 
fore, d x G xb must be of second order in d x at most, which 
implies that G xb involve, at most, first-order derivatives 
with respect to x. Of course, the reasoning runs as 
well for the y— derivatives, the z— derivatives, and the 
t— derivatives, so that G ay has not second y— derivatives, 
G az has no second z— derivatives and G th have no sec- 
ond time derivatives. In fact, traditionally this reasoning 
has been implicitly used to argue that the constraints 
G tb — are the only conditions needed at the initial slice 
to specify a solution to the Einstein equations via the 
initial- value problem [bH ]. 

We conclude that G a be b = g a cG cx = contain as much 
knowledge about the solution of the Einstein equations 
in the interior as the boundary values may be expected 



to "know". For this reason, we think that imposing the 
boundary conditions © on the incoming characteristic 
fields is necessary in order to evolve a solution to the evo- 
lution equations from constrained initial data that will 
satisfy the four constraints at a later time slice as well. 

To illustrate the point, we now consider the Einstein- 
Christoffel formulation of the 3+1 equations By 
defining first-order variables as the following 18 linearly 
independent combinations of the space-derivatives of the 
metric: 

fkij = r(ij)fc + Ikil m F[lj] m + 7fcj7 im r \li] m (8) 

where T k ij are the Christoffel symbols of 7y, and choos- 
ing the lapse function as a = Q^/7 with Q assumed arbi- 
trarily prescribed a priori, the 3+1 equations can be put 
in the following equivalent form [l6| : 



7y = ~2aK i3 (9a) 

K tj + al kl dif kl] = a{~/ kl (K k iK i3 - 2K ki K l:j ) + "f kl "f mn ^f kml f[ln]j + ±fkm[nfl]ij ~ fikmfjln 
+8f(ij)kf[ln]m + 4:fkm(ifj)ln ~ 8fklifmnj + 20/fe;( i /j) mn — Hfiklfjmn) 

-did, hxQ-di In Qdj In Q + 2 1%n kl 1 mn {f kmn d l In Q - f kml d n In Q) 
+7 fcZ [(2/ fo)fe - fkii)dl lnQ + ±fkl{idj) InQ - 3(/j«fy InQ + f jH di InQ)]} (9b) 
fkij + ad k K tJ = a{7 mn [4if fe( j/ J ) mn - 4f mn (iKj)k + K lJ (2f mnk - 3f kmn )} 

+ 2'y mn Y q [Kmp{l k ( l fj) qn - ^Iqn{ilj)k) + lk(i K j) m (8f n pq ~ Gfpqn) 

+K mn (Af pq (^ j)k - 5j k(i f j)pq )} - KijdklnQ 

+2 1 mn (K m{llj)k d n In Q - K mnlk{l d 3) In Q)} (9c) 



with the constraints 



C = 



c, 



1 tl 1 kl {2{d k f lJ i - difjki) + K ik Kji - K l3 K H 

run 

+fijk (Sfrnln - 20f lmn ] } = (10a) 
7 {l mn [Kik{%flmn - 2fmnl) ~ K km f a . n ] 

+d t K kl - d k K a ) = (10b) 

dklij - 2 fkij + ^l lm {flm{ilj)k - 7k(ifj)lm) = 

(10c) 



to be imposed only on the initial data. Here C and Ci 
are the Hamiltonian and momentum constraints, respec- 
tively, and Cij k represent the definition of the additional 
18 first-order variables f k ij needed to reduce the ADM 
equations to full first-order form (they represent JSJ) in- 
verted for 7y,fe in terms of fkij)- No other constraints are 
needed to pick a solution of the Einstein equations out 
of the larger set of solutions of the Einstcin-Christoffcl 
equations. 

In this language, the boundary conditions (j6bl6d(l re- 
quire the derivatives of the lapse a and the components 
R x , R x and R x — (1/2)R expressed in terms of fkij, the 



point being that there is a way to express their occur- 
rence so that no a;— derivatives of any f k ij appears. As 
they stand, they remain as evolution equations for the 
components K X ,K X and Kv + K* respectively. However, 
(|6a|l takes the following form 



G tb e b 



rx 
J 7 



jy-rn j'X 



f m + a{2K l (4/ m 
Kd x a + K xl dia} 



3f 



i J 
(11) 



where f k j = dt(^ km fmWy )■ This is an evolution equa- 
tion for a particular linear combination of the compo- 
nents Of fkij. 

The next step is to figure out which variables these 
four equations should prescribe at the boundary. We 
start by recalling the characteristic fields of the problem 
with respect to the unit vector = j xi / ^f^ xx which is 
normal to the boundary x — xq for the region x < xq. In 
this case, we have that 



+ TJ3 



U =K 



(12) 



are six characteristic fields that are outgoing at the 
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boundary, and 



-m =k 



ex j 



(13) 



are another six characteristic fields that are incoming at 
the boundary. The remaining characteristic fields are the 
six components of f y i J the six components of f z ^ and 
the six components of 7^, all 18 of which travel upwards 
along the boundary (they have vanishing characteristic 
speeds). 

Inverting to find the fundamental fields in terms of the 
characteristic fields we have 



K 



/V = 



-W 



(14) 
(15) 



One can immediately see that the boundary conditions 
(|6b|l and ({6*cjl turn into evolution equations for ~Uy and 
~Ug assuming the outgoing fields + Uy and + U Z are 
known, respectively. In fact, these equations yield the 
values of ~Uy and linearly in terms of the outgoing 
fields and in terms of sources known from the previous 
time step. So the transverse components a = y, z of the 
projection G a be b = are quite welcome and are actually 
necessary and consistent with the initial-boundary value 
problem of the Einstein equations. 

Turning now to the longitudinal component a = x of 
G ab e b = 0, which is Eq. JMJ, we see that, because K — 
K'i — Kjj + K*, the equation takes the form 



On the other hand, because 



= 0. (16) 



r 



/V + /"": 



/v - r 



(17) 



then the time component a = t of G a fe e = 0j which is 
Eq. Qllfl. takes the form 



{ + uy+ + ui)-{-uy+-u* z ) 



_ fy x 
J V 



■/%*+... = 0. (18) 



Both equations thus involve the same incoming and out- 
going characteristic fields: - Z7* + ~U% and + Uy + + U* . 
This means that they determine both the outgoing and in- 
coming fields + U-j+ + U z and ~U^+~U*. One way (by no 
means the only one) to make this point explicit is the fol- 
lowing. Subtracting the two equations with appropriate 
factors one obtains an evolution equation for ~U% + ~U* , 
which is welcome as a prescription of boundary values 
that are consistent with the Einstein equations: 



2aG x be 



;G tb e b = 



-C/* + -E£ + ... = 0, (19) 



which also involves time derivatives of some of the fields 
that travel along the boundary, but this is fine. In fact, 
we have shown already that at least in the case of 



the Einstein-Christoffel formulation with spherical sym- 
metry, this resulting boundary evolution equation is be- 
ing used in numerical integration and is referred to as a 
constraint-preserving boundary condition 9] . 

So far, the projection of the Einstein equations in the 
direction normal to the boundary is seen as providing 
prescriptions for the boundary values that are necessary 
in the light of the initial-boundary value problem of this 
strongly hyperbolic formulation. Any solution to the 
ten Einstein equations must satisfy such boundary equa- 
tions. Therefore, failure to impose such conditions on 
the boundary would result in a solution to the evolution 
equations with constrained initial data which would not 
solve all four of the constraints at a later time slice. 

There remains the question of the fourth boundary 
equation implied by adding Gtbe b = and G x \,e = 
with appropriate factors in order to eliminate the occur- 
rence of the incoming fields. This equation is linearly 
independent of the other three, therefore it will not be 
satisfied as a consequence of imposing the other three. 
One can see by inspection that this is an evolution equa- 
tion for the outgoing field + U y + + U*: 



2aG x be 



= 0, (20) 



which involves also the time derivatives of some of the 
fields that travel along the boundary, but it does not 
involve the time derivative of any incoming characteris- 
tic fields. This equation is not "necessary" in the sense 
that the outgoing fields are already determined by prop- 
agation from their initial values. However, this equation 
must be consistent with the outgoing fields determined 
by the initial values, otherwise one would not be solving 
the full set of Einstein equations. Therefore, it is actu- 
ally necessary for this equation to be satisfied. In fact, in 
principle, one may simply use this equation for those out- 
going characteristic fields involved, with the confidence 
that this boundary prescription is consistent with the 
initial values, so long as the initial values satisfy the con- 
straints with sufficient accuracy. The actual numerical 
implementation of this boundary condition is a subject 
for subsequent study and lies beyond our current scope. 

In summary, in the case of the Einstein-Christoffel for- 
mulation, the four "Einstein boundary conditions" de- 
fined as G a be b — up to linear combinations constrain 
the values of three of the six incoming characteristic 
fields. The other three incoming characteristic fields are 
free, since there are no other available conditions on the 
boundary. In fact, the 18 first-order constraints Cijk - 
which are needed in order to pick an actual solution of the 
Einstein equations out of the larger set of solutions of the 
Einstein-Christoffel equations- are preserved trivially by 
the evolution. This can be seen easily by taking a time 
derivative of C^k and using the evolution equations to 
eliminate the occurrence of time derivatives of the funda- 
mental variables in the right-hand side. The reader may 
verify that this leads to C^k — Fijk(C,Ci,Cijk), that 
is: ordinary differential equations for the propagation of 
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Cijk- This means that the 18 first-order constraints are 
trivially satisfied on the boundary by virtue of the ini- 
tial data and the evolution equations, and, consequently, 
they may not prescribe conditions on the incoming fields. 
There thus remain three incoming fields that can be spec- 
ified in any arbitrary fashion for a solution to the Einstein 
equations. They may, for instance, be "frozen" by set- 
ting their time derivatives to zero, as is done by some 
authors [ll . 

We have thus obtained a complete set of boundary con- 
ditions for the Einstein-Christoffel formulation from the 
projection of the Einstein tensor normally to the bound- 
ary and the fact that the 18 constraints that realize the 
first-order reduction are preserved by ordinary differen- 
tial equations. The incoming fields that are prescribed 
by the three boundary equations carry, in a sense, the 
vanishing of G a be b into the interior, where the evolution 
equations hold. Since G a be b are linear combinations of 
the evolution equations and the constraints, it follows 
that the constraints will be satisfied in the interior. Thus 
the boundary conditions preserve the constraints. 

The argument makes no use of the auxiliary system 
of propagation of the constraints implied by the evolu- 
tion equations. Still, the resulting boundary conditions 
appear to be completely consistent with constraint prop- 
agation in the following sense. In the three-dimensional 
case that interests us here, if the evolution equations are 
satisfied, the constraint functions propagate according to 
a strongly hyperbolic system [T^. In this auxiliary sys- 
tem, three characteristic constraints are incoming, three 
are outgoing and the rest are "static". In the context 
of this auxiliary system, the three incoming constraints 
must be set to zero along the boundary. One may thus 
infer that there will be three conditions on the boundary 
values of the fundamental fields that are necessary and 
sufficient to preserve the constraints. Since we have three 
Einstein boundary conditions and no other conditions on 
the boundary, the Einstein boundary conditions and the 
three conditions that arise from constraint propagation 
must be equivalent in some definite sense. We argue that 
the equivalence is through linear combinations of the evo- 
lution equations. In the spherically symmetric case we 
have shown that this is precisely the case ^l|. Work is 
currently in progress to demonstrate explicitly the equiv- 
alence of the two sets of conditions in the three dimen- 
sional case and will appear elsewhere. Interestingly, since 
constraint propagation is at the basis of Calabrese et 
al's "constraint-preserving" boundary conditions [l^, it 
would thus appear that projecting the Einstein equations 
normally to the boundary would shortcut the procedure 
leading to "constraint-preserving" boundary conditions, 
providing an exact nonlinear version of the conditions 
and a geometric basis for them at the same time. 

We emphasize that there is a very subtle issue of con- 
sistency between the constrained initial data traveling 
towards the boundary, and the conditions G a be b = on 
the boundary that apply to the outgoing fields. The is- 
sue turns critical in numerical simulations because of the 



inaccuracies involved in satisfying the initial constraints. 
Clearly, the outgoing fields propagating out from the ini- 
tial data will carry and perhaps increase the effects of the 
initial errors, and there is a good chance that they will fail 
to satisfy l|20|l. the odds of failure increasing with time. 
This issue may have large consequences for achieving long 
time numerical evolution. This observation applies espe- 
cially to numerical simulations implementing black hole 
excision 0|- As far as we know, no consistency issues in 
strongly hyperbolic systems with constraints have been 
addressed in the reference literature on finite-difference 
methods such as [l3|. 

Interestingly, the other constrained system of equa- 
tions in physics which is oftentimes used as a model for 
the Einstein equations -that is, the Maxwell equations- 
does not have this issue. When the Maxwell equations 
d a F a i, = Jb are written in terms of a vector potential 
F a b = d[ a A b ] , the other half of the system d[ a F bc ] = is 
identically satisfied and the time component Aq = <p of 
the vector potential may be considered arbitrary (just 
as the lapse function is arbitrary insofar as the 3+1 
equations are concerned). It is a straightforward ex- 
ercise to put the equations in "ADM form" by defin- 
ing a variable Ki = A4, and to show that the projec- 
tion e b d a d[ a A b ] — e b Jb for e b = S b turns out to be ex- 
actly one of the evolution equations. Exactly the same 
thing happens in the Lorenz gauge. The underlying rea- 
son being, of course, that there is only one constraint 
C = n b d a d[ a A b ] — n b Jb with n b = <5(j which is preserved 
by the evolution by means of the trivial equation C = 0. 

A separate issue not treated here is whether or not the 
"Einstein boundary conditions", that is: G a be b = up to 
linear combinations, would lead to a well-posed initial- 
boundary value problem for the Einstein-Christoffel for- 
mulation in the sense that the boundary conditions pre- 
serve the well-posedness of the initial value problem. In 
this respect, Calabrese et al ^8| show how to obtain a 
subset of the constraint-preserving boundary conditions 
that lead to an overall well-posed initial-boundary value 
problem -including constraint propagation- for the lin- 
earized "generalized" Einstein-Christoffel system. The 
main difference between the Einstein-Christoffel system 
and its modified counterpart (aside from the fact that 
the variable fkij is defined differently) is that the con- 
straints of the modified system propagate in a symmet- 
ric hyperbolic manner, a difference that does not affect 
the well-posed character, so one may expect that a sub- 
set of well-posed boundary conditions can be found in 
the Einstein-Christoffel case as well. Furthermore, the 
"Einstein boundary conditions" clearly exist and can be 
written down for any formulation of the Einstein equa- 
tions, including those that are symmetric hyperbolic and 
that have symmetric hyperbolic constraint propagation, 
by simply expressing Eqs. © in terms of the fundamental 
variables of the problem. There are preliminary indica- 
tions that, in such cases, linear combinations of G a be b = 
exist which are well posed in the sense of [l8|. 

But most of the numerical relativity effort is cur- 
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rently driven by formulations that are not strongly hy- 
perbolic (see, for instance, [H EJ |H |H |M H |2(| ) in 

which case sensible boundary conditions are also needed, 
whether or not they are well-posed. For those formula- 
tions, the constraint-preserving boundary conditions of 
Calabrese et al may not even be defined, yet the pro- 
jections G a bS b = evaluated on the boundary are al- 
ways available, are easy to write in exact non-linear form, 
and provide physical boundary conditions that are clearly 
necessary in order to enforce the constraints, the question 
of well-posedness becoming irrelevant. 

There is the question of why should numerical simu- 
lations use these boundary equations at all if they can 
do just fine by simply "pushing the boundaries out" - 
even though doing so increases the computational re- 
quirements considerably. In this respect, we point out 
the existence of the following conjecture based on the 
fact, observed by many authors, that the run time of nu- 



merical simulations is noticeably shortened by growing 
constraint violations. The conjecture, yet to be proven 
rigorously, is that the run time might be extended by 
efficiently controlling the constraint violations. Since we 
have shown that the boundaries are intimately connected 
with constraint propagation, if the conjecture turns out 
to be true then the "Einstein boundary conditions" may 
have a role to play in extending the run time of numerical 
simulations. 
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